Density Functional Study of Ternary Topological Insulator Thin Films 
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Using an ah initio density functional theory based electronic structure method with a semi-local 
density approximation, we study thin-film electronic properties of two topological insulators based 
on ternary compounds of Tl (Thallium) and Bi (Bismuth). We consider TlBiX2 (X=Se, Te) and 
Bi2X2Y (X,Y= Se,Te) compounds which provide better Dirac cones, compared to the model binary 
compounds Bi2X3 (X=Se, Te). With this property in combination with a structurally perfect bulk 
crystal, the latter ternary compound has been found to have improved surface electronic transport 
in recent experiments. In this article, we discuss the nature of surface states, their locations in 
the Brillouin zone and their interactions within the bulk region. Our calculations suggest a critical 
thin film thickness to maintain the Dirac cone which is significantly smaller than that in binary 
Bi-based compounds. Atomic relaxations or rearrangements are found to affect the Dirac cone in 
some of these compounds. And with the help of layer-projected surface charge densities, we discuss 
the penetration depth of the surface states into the bulk region. The electronic spectrum of these 
ternary compounds agrees very well with the available experimental results. 

PACS numbers: 71.15.Dx, 71.18.+y, 73.20.At, 73.61.Le 



I. INTRODUCTION 



Three dimensional (3D) Topological band Insulators 
(TI) have attracted considerable attention from the con- 
densed matter and device physics communities for the 
novel electronic surface states they support and the host 
of unusual responses to external fields 1 -. The basic un- 
derstanding of their electronic structure, using Density 
Functional Theory (DFT) is a continuing quest in the 
literature. These studies have contributed to the growth 
of the field and helped with interpreting experiments. In 
recent years, binary 3D TI materials Bi2Se3 and Bi2Te3 
have emerged as model systems for numerous experi- 
ments which focused on exfoliation and the molecular 
beam epitaxy growth of thin films and for probing the 
electronic structure by optical techniques^ and device 
fabrications-. The theoretical works, based on density 
functional theory or tight-binding calculations, focused 
on elucidating the structure-property relations, predict- 
ing new TIs, and mapping the band structures 4 . One 
of the concerns for the use of binary TI materials in de- 
vices is the intrinsic n(or p)-type vacancies in the bulk 
crystals of Bi2Se3 (Bi2Tes) which make the conduction 
through bulk states dominate transport experiments 5 
and the hexagonal warping effect that makes the Dirac 
cone anisotropic, especially in the conduction band re- 
gion. This anisotropic nature of the Dirac cone leads 
to interesting spin-textures (or spin-momentum locking) 
of surface states with possible applications in spintron- 
ics and quantum information processing^. Recently, two 
promising 3D TI materials, TlBiX 2 (X= Se,Te) and 
Bi2X2Y (X,Y = Se,Te) have been predicted to have near 
perfect Dirac cones, in terms of less entanglement of bulk 
and surface states, and experiments have supported these 
predictions.— ~—. Moreover, Bi-based ternary compounds 
offer high bulk resistivity, due to the structurally perfect 
nature of the crystals, so surface transport is enhanced. 
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FIG. 1: (Color online) Schematic of the hexagonal bulk struc- 
tures of representative ternary compounds (a) TlBiSe2 and 
(b) bulk Bi2Se2Te derived from their corresponding bulk trig- 
onal structures whose three primitive vectors are denoted by 
ti, t2 and t3. The trigonal structure is part of the larger 
hexagonal structure shown by dashed lines. In both the com- 
pounds, the hexagonal cell contains three times the number 
of atoms compared to that in the trigonal structure. In (b), 
the five atomic layers, forming a quintuple layer are shown 
in the shaded square region, (c) The first Brillouin Zone of 
the bulk trigonal structure with four time-reversal invariant 
points r, Z, F and L are shown. 



In view of these promising advances, a comprehen- 
sive theoretical study of the thin-film structures of these 
ternary 3D TI is necessary to help better understand 
these materials and to help design future experiments 



2 



to verify these effects. These materials can have intrinsic 
size limits which can protect the metallic nature of the 
surface bands, the surface states can spread inside the 
bulk region, and atomic rearrangements in thin layers can 
have profound effect on the Dirac cone itself. We address 
these issues in this paper using a DFT-based electronic 
structure method and compare our results with available 
experimental results as well as with the results of binary 
Bi-based TIs. Our studies will have implications for the 
understanding of topological surface states in ternary TIs 
and their intended applications. 

Our paper is organized as follows. In Section II, we de- 
scribe the bulk crystal structures of TlBiX 2 (X=Se,Te) 
and Bi 2 X 2 Y (X,Y=Se, Te), their conduction band (CB) 
and valence band (VB) structures and the computational 
method used for this study. In Section III we present 
the thin film electronic structure of these materials and 
discuss the role of atomic relaxations or rearragemcnts, 
resulting from thin film formation from the bulk crys- 
tal, on the shape and size of the Dirac cone in the bulk, 
and predict critical film thicknesses required to maintain 
the Dirac cone and degree of surface state extension into 
the bulk region which have yet to be determined experi- 
mentally for these class of materials. We compare these 
results with those available for binary TI compounds. 
Finally, in Setion IV, we present our summary and con- 
clusions. 



II. COMPUTATIONAL METHOD AND BULK 
BAND STRUCTURES 

This section details the computational method used, 
the choice of computational parameters, the bulk crys- 
tal structures and the resulting electronic band struc- 
ture. The bulk crystal structures of both Tl and Bi- 
based ternary compounds are similar to the binary com- 
pounds Bi 2 X3(X=Se,Te). Both compounds have a trigo- 
nal structure with covalently bonded alternating cations 
and anion layers stacked along the crystallographic z- 
direction. However, there is one important difference as 
compared to the binary compounds: The atomic layers of 
these Tl-based structures are not arranged in a quintuple- 
like (QL) order, suggesting covalent bonding between the 
unit-cells, and therefore, that exfoliation techniques, used 
to peel binary TI flakes from the corresponding bulk crys- 
tals, cannot be used for extractng Tl-based TI thin films. 
The layers are arranged in the order Tl-Se(Te)-Bi-Se(Te). 
The trigonal unit cell contains four atoms (as opposed to 
five in binary Bi-based TIs) with the lattice parameter a 
= 0.7887 nm(0.8263 nm) and the angle between the lat- 
tice vectors spanning the lattice of a~31.4°(~31.8°) for 
TlBiSc 2 (TlBiTe 2 The hexagonal cell, formed from 
the trigonal, consists of twelve atomic layers with the 
number of atoms tripled (Fig. 1(a)). The lattice param- 
eters of the hexagonal cell are a=0.4264 nm(0.4534 nm) 
and c=2.2478 nm(2.35I2 nm) for TlBiSe 2 (TlBiTe 2 ) and 
the layer sequence in this case is same as the trigonal 



cell. Compared to the corresponding binary compound, 
the Bi-based ternary compound Bi 2 Se 2 Te (Bi 2 Te 2 Se) has 
one of the Se (Te) layers replaced by a Te (Se) layer, and 
has the layers arranged in the order Te(Se)-Se(Te)-Bi- 
Bi-Se(Te). The QL-like layer arrangement in these com- 
pounds facilitates exfoliation-like methods which can be 
used to peel thin films from the bulk- crystals. The lat- 
tice parameters of the bulk trigonal cell is almost same 
for Bi 2 Se 2 Te (Bi 2 Te 2 Se): a = 1.0046 nm(1.0255 nm) and 
a ~ 24.2° (~ 24.1°) 2 i The hexagonal cell, built from the 
trigonal structure, contains fifteen atomic layers with the 
lattice parameters a = 0.422 nm(0.428 nm) and c = 2.92 
nm(2.99 nm) for Bi 2 Se 2 Te (Bi 2 Te 2 Se) (Fig. 1(b)). 

We used a DFT-based electronic structure method 
with projector-augmented wave basis^ and a general- 
ized gradient approximation to the exchange-correlation 
potential^ for computing the electronic properties of the 
bulk and thin films of both of the compounds. Spin-orbit 
coupling (SOC) was invoked in the calculation as imple- 
mented in the numerical method^. For the bulk struc- 
tural optimizations of both the compounds, a kinetic en- 
ergy cut-off of 400 eV and a non-orthogonal k-point mesh 
along the reciprocal lattice vectors of 9 x 9 x 9 in the first 
Brillouin zone (BZ) were chosen. Since previous struc- 
tural studies of bulk trigonal structures reported agree- 
ment with experimental lattice parameters^, we chose to 
relax only the internal parameters, namely the atomic po- 
sitions, keeping the lattice constants fixed to experimen- 
tal values. The total energy is assumed to be converged 
when all the components of Hellman-Feynman forces on 
each ion are smaller than the threshold 0.001 eV/A. Con- 
vergence of computed properties were carefully checked 
with respect to large energy cut-off, k-point mesh and 
larger force threshold. For both the compounds, the com- 
puted interlayer distances are quite close to the experi- 
mental worksi^. These computed distances suggest van 
der Waals type of bonding between the adjacent QLs in 
ternary Bi-based TI but not in Tl-based TIs. With these 
informations, we built the bulk hexagonal cell and the 
corresponding thin-films of both the compounds as de- 
tailed in the next section. 

The bulk band structures of TlBiX 2 (X= Se,Te), with- 
out and with SOC, agree well with other theoretical 
calculation o 15 ' 17 . However, with SOC, other calcula- 
tions predict TlBiTe 2 to be an indirect semiconductor 
with the CB minimum (CBM) at the T point and the 
VB maximum (VBM) on the line joining the L and Z- 
points. The high-symmetry points are labeled in Fig. 
1(c). Recent experimental studies argue TlBiTe 2 to 
be a semimetal with an negative energy gap of 20 
meV— . Our results agree well with this experimental 
work (Fig. 2(a)). TlBiSe 2 has a direct gap at F of 124 
meV(235 meV) without(with) SOC. Figure 2(b) and (c) 
show the bulk band structures of Bi 2 X 2 Y (X,Y=(Te,Se) 
or,(Se,Te)) with SOC, which agree quite well with a re- 
cent DFT calculation^. The CBM is at the T-point and 
the VBM is between Z and F-points for both the Bi- 
based TIs. The computed indirect gap are 157 meV and 
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FIG. 2: (Color online) Band structure of bulk trigonal struc- 
ture of (a) TlBiTe 2 and (b) Bi 2 Se 2 Te and (c) Bi 2 Te 2 Se with 
spin-orbit coupling, shown along high symmetry dirctions in 
the bulk BZ. We obtained a negative indirect gap in TlBiTe 2 
and indirect gaps in Bi-based ternary TIs. Energy states 
mainly from p z orbitals of Bi, Se and Te are marked on the 
top of bulk bands and the procedure to obtain these surface 
bands from crystal wave-functions is discussed in the text. 



272 meV, for Bi 2 Se 2 Te and Bi 2 Te 2 Se, respectively. We 
note that the gap is direct at the T-point without includ- 
ing SOC with the gap values of 633 meV and 76 meV, 
for Bi 2 Se 2 Te and Bi 2 Te 2 Se, respectively. 

In the literature, emergence of these topological sur- 
face states within the bulk band gap of binary 3D TIs 




FIG. 3: (Color online) Band structure of bulk Bi 2 Te 2 Se along 
high symmetry directions (a) without spin-orbit coupling, and 
(b) with spin-orbit coupling. Energy states formed mainly 
from p z orbitals of Bi and Te are marked on the top of bulk 
bands and band inversion processes are clearly shown. 



was associated with the bulk band inversion when SOC 
is switched oi>k We tested this conjecture by plotting 
the specific orbital contributions from a given atom as a 
function of k and band on the band structure diagrams 
of Figs. 2 and 3. These specific contributions are nor- 
malized with respect to the contributions from all the 
orbitals in both the bulk Tl- and Bi-based ternary com- 
pounds. Since the bulk CBM and VBM at T consist 
mainly of p z orbital*^ of the Bi and Se(Te) atoms, re- 
spectively for Tl-based compounds, this state was cho- 
sen for the orbital projection studies. In Tl-based TIs, 
a cut-off contribution percentage is chosen which gives a 
reasonable picture of band inversion in these figures. If 
the calculated p z orbital contribution of a specific atom 
at a certain state is greater than a given cut-off percent- 
age, we considered this state as mainly originating from 
that atom. This cut-off was set at 50% and 25% for 
Se(Te) and Bi atoms, respectively. A slight change in a 
cut-off percentage, below or above these choices, resulted 
in an insignificant change for the CBM and VBM bands 
at the T-point, which is region of interest to study the 
band inversion process. For example, the p z orbital con- 
tributions on the band structure of TlBiTe 2 in Fig. 2(a) 
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FIG. 4: (Color online) Schematic diagram of thin film struc- 
tures of (a) TlBiSe2 obtained by stacking 39 layers along the 
z-direction. The atomic layers close to surfaces are in the 
order of Se-Bi-Se-Tl, on both sides of the film, (b) Two- 
dimensional Brillouin zone of the (111) surface of the bulk 
TlBiSe2 with three time-reversal invariant points T, M, and 
K 



remain similar near the CBM and VBM around the T- 
point with the cut-off percentage of 45% and 22.5% for 
Te and Bi atoms, respectively (Figures not shown). A 
similar band inversion effect was also seen in Bi-based 
TIs (Bi 2 Te 2 Se and Bi 2 Se 2 Te). This effect is shown only 
for one of the compound Bi 2 Te 2 Se (Fig. 3(a) and (b)). 
The other compounds show similar effects (Figures not 
shown) . We gave an equal weight of 30% to both Bi and 
Te(Se) orbital contributions since the two outermost lay- 
ers of Se(Te) and Bi, out of the total of five atomic layers 
within each QL, mainly particpate in forming the VBM 
and CBM. 



III. THIN FILM BAND STRUCTURES AND 
SURFACE STATES 

In this section we discuss the surface states of ternary 
TIs. First we discuss the interesting nature of surface 
states in Tl- and Bi-based TIs. We compare our results 
with binary Bi-based TIs, wherever possible. We choose 
kinetic energy cut-off of 400 eV and k-mesh size of 9 x 
9 x 1 on the surface BZ of the hexagonal supercell for 
computing surface band structure of both Tl- and Bi- 



based ternary TIs. 



A. Surface states in Thallium-based ternary 
compounds 

The thin films of Tl-based TIs are built from the bulk 
hexagonal structure with the atomic layers stacked along 
the crystallographic z direction and a net vacuum size of 3 
nm above the top layer and below the bottom layer within 
a periodic simulation region. Recent theoretical^ stud- 
ies suggest contributions to the surface band structure 
from the dangling bonds arising from the surface termi- 
nations, and that these dangling states appear along with 
the topological surface states. Among four surface ter- 
minations, with the atomic layer ordering as Tl-Se(Te)- 
Bi-Se(Te), Se(Te)-Bi-Se(Te)-Tl, Bi-Se(Te)-Tl-Se(Te) and 
Se(Te)-Tl-Se(Te)-Bi, the one with Se(Te) beneath the Tl 
layer is argued to have minimal surface dangling bond 
densitys-. Therefore we chose this termination in all fur- 
ther thin film studies. We note here that experiments 
do not observe any signature of dangling bond states in 
the optical spectral. And, at least for this termination, 
we do not find any signature of the dangling bond state 
either, at least, within the Dirac cone region. Different 
thicknesses of the thin films of TlBiSc 2 , ranging from 
23 to 39 layers, corresponding to thicknesses of 4 to 7 
nm and 23 to 31 layers for TlBiTe 2 with the maximum 
thickness value of 5.8 nm were considered. Figure 4(a) 
shows a representative supercell structure of 39 layers of 
TlBiSe 2 . We consider the relaxation of atomic positions 
in the z-direction since it is argued that the positions of 
Se(Te) atoms are critical in determining the topological 
nature of the surface states^. The positions are assumed 
to be optimized when the ^-component of the forces are 
smaller than the threshold value 0.015 eV/A. 

As for the Bi-based binary compounds, Tl-based 
ternary compounds also show thickness dependent elec- 
tronic structure. The Dirac cone is preserved in TlBiSc 2 
with 39 atomic layers, corresponding to a thickness of 
7 nm (Fig. 5(a)), and in TlBiTe 2 with 31 layers of a 
thickness 5.8 nm (Fig. 5(b)), where red circles denote 
surface state contributions and the solid black lines rep- 
resent the bulk bands. We reproduce the interesting ex- 
perimental features of TlBiSc 2 7 and TlBi 2 Te±i. The bulk 
VB (BVB), in TlBi 2 Se has two maxima between T and 
M points in the hexagonal BZ (Fig. 4(b)) and the bulk 
CB (BCB) has a minimum at M. The Dirac cone is found 
to exist inside the bulk band gap, at the F-point of the 
BZ. These features are clearly seen in our computed band 
structure (Fig. 5(a)). Our computed band structure of 
TlBiTe 2 shows that the Dirac point is below the BVBM 
(Fig. 5(b)). The BCBM is placed at f and the BVBM 
along the line joining T and M, producing an indirect 
gap of 83 meV. Another BVBM emerges along the line 
joining F to K. These features including the indirect gap 
size, positions of BVB, BCB and the Dirac point are con- 
sistent with the recent experimental studyii. 




FIG. 5: (Color online) Band structures of (a) TlBiSe2 film 
with a thickness of ~7.0 nm (39 layers) and (b) TlBiTe2 film 
with a different thickness of ~5.8 nm (31 layers). Thinner 
films (23, 27, 31, 35 layers for TlBiSe2 and 23, 27 layers for 
TlBiTe2) show nonzero band gaps at F (Figures not shown). 
Orbital contributions from first few layers in each thicknesses 
are marked with red circles. 



For thin film thicknesses corresponding to fewer than 
39 (31) atomic layers in TlBiSe 2 (TlBiTe 2 ), our study 
suggests that a band gap opens in the otherwise metallic 
surface states, much like those in binary Bi-based TIs 19 . 
With decreasing thicknesses, the size of the gap seems 
to increase indicating the increased interactions between 
surfaces states that exist on two opposite surfaces of the 
TI (Figures not shown). The gap values are summarized 
in Table I. There are no experimental reports, so far, 
on the critical film thickness necessary to preserve the 
metallic bands for Tl-based ternary TIs which our theo- 
retical study predicts. It is interesting to note that the 
predicted induced gaps in Tl-based TIs, for the same film 
thickness, are larger than the binary Bi-based TIs22,. 

To understand the origin of the gap at the Dirac point 



0.4 



0.3 



> 0.2 
CD 

w 



-0.1 



-0.2 



FIG. 6: (Color online) Band structure of TlBiSe2 film with a 
thickness of ~5.6 nm (31 layers). A nonzero band gap exists 
at f due to the coupling between the top surface CB and the 
bottom surface VB and between the top surface VB and the 
bottom surface CB. No splitting is observed away from the 
surface state band edges. 
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FIG. 7: (Color online) Layer-projected charge density, con- 
tributed by surface states in the neighborhood of the F-point 
for (a) TlBiSe2 with a thickness of ~ 7 nm, corresponding to 
39 atomic layers and (b) TlBiTe2 with a thickness of ~ 5.8 
nm, corresponding to 31 layers. The color labeling scheme is 
consistent with that in Figure 1 (a). 
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TABLE I: The induced band gap at the time-reversal invari- 
ant point r AEf (in eV) for thin-films with various thick- 
nesses. 



Number of Layers 


23 


27 


31 


35 


39 


AEf 
TlBiSe 2 


0.059 


0.042 


0.023 


0.015 


0.000 


AEp 
TlBiTe 2 


0.018 


0.009 


0.000 


NA 


NA 




1 4 8 12 16 20 24 27 

Layer 



FIG. 8: (Color online) Layer-projected relative charge density 
for thinner samples of (a) TlBiSe2 with thickness of ~ 5.6 nm 
corresponding to 31 layers and (b) TlBiTe2 with thickness of 
~ 5 nm corresponding to 27 layers. The color labelling scheme 
is consistent with that in Figure 1(a). 



with decreasing film thickness, we first map out the sur- 
face state contributions in the band structure from atoms 
within first few layers of the film, and from all the or- 
bitals, relative to the total contributions of all orbitals 
from all the layers of the thin film. From the top or bot- 
tom surface of the film, 8(6) atomic layers were chosen for 
TlBiSc2(TlBiTe2). A larger number of atomic layers are 
used for TlBiSe2 than for TlBiTe2 because a larger num- 
ber of atomic layers is required to preserve the metallic 
band structure because the surface states are less local- 
ized to the surface region in TlBiSe2- These choices also 
helped us to choose the same cut-off percentage for both 
the Tl-based compounds in estimating the surface state 
contributions. For film thicknesses corresponding to 23, 
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FIG. 9: (Color online) Schematic diagram of (a) Bi2Se2Te 
and (b) Bi2Te2Se thin film structures obtained by stacking 
4QLs along ^-direction. 



27, 31, 35 and 39 atomic layers in TlBiSe2, respectively, 
these cut-offs were varied from 70% - 50%, in steps of 5% 
decrease. For TlBiTe2, 23, 27 and 31 atomic layers were 
considered with cut-offs of 70%, 65% and 60%. Examples 
results of these orbital contributions are marked on the 
top of the total band structure in Fig. 5(a) and (b). 

After estimating the distribution of surface states in 
the thin film structure, we computed the valence charge 
density resulting from these surface state wave-functions 
in both the compounds. We focused on surface states 
around the Dirac-point, in an energy window chosen such 
that same number of wave-functions are used for build- 
ing the surface charge density for most relevant thick- 
nesses considered in this study. The wave-functions at 
three closely spaced k-points including the T-point were 
chosen. The charge density is computed layer-wise. We 
define the relative charge density, to be the ratio of peak 
charge density in a particular layer to the peak charge 
density in all layers. This quantity is plotted in Figs. 
7(a) and (b). Our studies suggest exponential decay of 
the charge density in the bulk region but with slower de- 
cay in TlBiSe 2 than in TlBiTe 2 . In about 12-14 atomic 
layers, from top or bottom surfaces of the film (corre- 
sponding to 2.3 nm), the surface state contributions de- 
cay to almost zero in TlBiTe2, whereas, the surface states 
extends up to 16-18 atomic layers (~3 nm) in TlBiSc2. 
With decreasing film thickness, the charge density accu- 
mulates in the middle of bulk region in both compounds 
(Figs. 8(a) and (b)) suggesting increasing interactions of 
surface states by overlap of their wave- functions inducing 
a gap at the Dirac point. However, while a band gap is 
opened in the surface state band structure, away from 
the surface state band edges band structure is largely 
unaffected as seen in Fig. 6 -i.e., it does not split into 
two bands as one would expect from resonant coupling 
between opposite surface states-suggesting that the split- 
ting is due to interband coupling between the top surface 
CB and the bottom surface VB and between the top sur- 
face VB and the bottom surface CB. 
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FIG. 10: (Color online) Band structures of Bi2Se2Te films 
with thickness of 3.6 nm corresponding to 4 QLs (a) without 
atomic relaxation and (b) with atomic relaxation. The Dirac 
cone is affected by relaxing the atoms. 



FIG. 11: (Color online) Band structures of Bi2Te2Se films 
with thickness of ~ 3.7 nm corrresponding to 4 QLs (a) with- 
out atomic relaxation and (b) with atomic relaxation. The 
Dirac cone is not affected by relaxation. 



B. Surface states of Bismuth-based ternary 
compounds 

We discuss the surface state properties of Bi2Se2Te and 
Bi2Te2Se in this section. The building blocks of these 
compounds are QLs, arranged in the order -Se(Te)-Bi- 
Te(Se)-Bi-Se(Te)-, and stacked along the crystallographic 
z-direction (Figs. 9(a) and (b)). Different film thick- 
nesses corresponding the different number of QLs were 
considered with a maximum of 4 QLs. The choice thick- 
ness is guided by the necessity to maintain the metallic 
nature of the surface bands. The computational param- 
eters used in this study are same as those for Tl-based 
TIs. We first construct the thin film from experimen- 
tal parameters^ and then study the thin-film properties 
without and with relaxing the atomic positions with the 
force convergence criteria of 0.05 eV/A. The procedure 
for extracting surface states from crystal wave-functions 



and the construction of layer- wise surface charge density 
is same as that for Tl-based TIs discussed in the previous 
section. 

The thin-film structure of Bi2Se2Te, without internal 
position relaxation, shows a more symmetric Dirac cone 
with metallic bands, as compared to the Tl-based TIs 
for thickness corresponding to 4 QLs (Fig. 10 (a)). The 
bulk band gap (of 313 meV) supports these novel states, 
and the Dirac point is well within this gap. Recently, 
a report of 3 QLs, as the minimum thickness, required 
to preserve the Dirac cone has appeared^. We find 4 
QLs to be necessary, which we also illustrate with the 
help of surface charge density calculations. However, in 
the absence of any experimental work, it is likely that 
the minimum thickness, carefully predicted by different 
theoretical methods, may fall in the range of 3QLs-4QLs. 
For film thicknesses below 4 QLs, a finite gap is induced 
at the Dirac point whose size increases with decreasing 
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(a) 




Layer 

FIG. 12: (Color online) (a) Layer-projected charge density 
for Bi2Se2Te without the atomic relaxation, contributed by 
surface states in the neighborhood of the F-point for thickness 
of ~3.6 nm corresponding to 4 QLs. Charge density resulting 
from (b) top and (c) bottom surface wave functions and (d) 
Bi2Fe2Se with the atomic relaxation for the thickness of ~3.7 
nm corresponding to 4 QLs. Although a charge pile up in the 
middle is seen in (a), the separate wavefunction contributions 
in (b) and (c) clearly show that surface charge density strictly 
vanishes after 2 QLs from both the top and bottom surfaces 
of the film. 



FABLE II: The induced band gap at the time-reversal invari- 
ant point F AEp (in eV) for thin-films with various thick- 
nesses. 



Number of QLs 


2 


3 


4 


AEp 

Bi2Se2Te without relaxation 


0.128 


0.038 


0.000 


AEp 

Bi2Se2Te with relaxation 


0.086 


0.017 


0.000 


AE r 

Bi2Te2Se without relaxation 


0.039 


0.037 


0.000 


AE r 

Bi2Te2Se with relaxation 


0.027 


0.037 


0.000 



the film thickness (Figures not shown). The size of the 
gap values are listed in Table II. With atomic relaxation, 
the Dirac cone is affected and its symmetry is lost (Fig. 
10(b)). The BVB lies along the line joining f and M 
in the hexagonal BZ, and its maximum is positioned al- 
most in line with the Dirac point. This hints at crucial 
role played by rearrangement of atomic positions either 
insitu in thin-film preparation procedure or in presence of 
external environments. In Bi2Te2Se, a significant part of 
surface state band structure falls below the BVBM and 
the Dirac point occurs in the energy neighborhood of the 
T-point. Atomic relaxation has almost negligible effect 
on the overall band structure (Figs. 11 (a) and (b)) in 
this case. Decreasing film thickness induces a finite gap 
(Table II). These studies suggest that the Dirac cone re- 
mains protected in thinner layers of ternary Bi-based TIs 
than in binary Bi-bsed TIs or than in Tl-bascd TIs. 

The origin of the critical thickness needed to main- 
tain the metallic nature of the surface states can be un- 
derstood by studying the layer-dependent surface charge 
densities. The procedure is same as that used for Tl- 
based TIs, except the charge density in the interstitial 
regions is also considered in order to locate the surface 
state spread in the bulk region as precisely as possible for 
thinner films thicknesses. For 4 QL thick Bi2Se2Te films, 
the combined charge density from the surface states of 
both sides of the film seems to be significant in the mid- 
dle of the bulk region (Fig. 12(a)). However, we sepa- 
rately computed the charge densities resulting from top 
and the bottom surface states. As seen from the plots 
(Figs. 12(b) and (c)), the density associated with surface 
states drops rapidly beyond 2 QLs, indicating negligible 
interaction between these opposite surface states in the 
4QL film, consistent with the observed protected Dirac 
cone. Relaxation of atomic positions is found to have no 
significant effect. In the 4QL films of Bi2 Te2 Se, even 
the combined contribution from both surfaces is small 
in our studies with and without atomic relaxation (Figs. 
12(d)). 
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IV. SUMMARY AND CONCLUSIONS 

We use a density functional based electronic structure 
method to study the thin him surface state properties of 
Tl and Bi-based ternary topologicial insulators. These 
studies predict that the Dirac cone remains protected in 
thinner layers of ternary Bi-based TIs than in binary Bi- 
based TIs or than in Tl-based TIs, which may be advan- 
tageous for certain applications. However, we also predict 
that in the atomically relaxed strucures, the Dirac cone 
of all but Bi2Se3 lie in or near the BVB. However, the 
Dirac cone of Bi2Se2Te remains outside of the BVB in 
the absence of relaxation, as might occur with dielectrics 
replacing the air gap. And, of course, the bonding to the 
dielectric could directly affect the surface states as well. 



Our computed results agree very well with experimen- 
tal results where available, while other predictions such 
as the critical him thicknesses required to maintain the 
Dirac cone, size of the induced gaps, and the extent of 
the spread of surface states in the bulk region point to 
further needed experimental work. 
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